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I ' Abstract 

r~| ! We report preliminary results from a study of the charmless exclusive semileptonic decay 

-K~^~^v based on the data collected at the T(4S') resonance using the BABAR detector at SLAC. 
, ■ The analysis uses events in which the signal B meson recoils against a B meson that has been 
\ reconstructed in a semileptonic decay B'^ D^*^^l~V. We extract the total branching fraction 
I B{B^ —>■ TT^i^v) = (1.03 ± 0.25stat. ± 0.13syst.) x 10^^ and the partial branching fractions in three 
bins of g^, the invariant mass squared of the lepton-neutrino system. From the partial branching 
fractions and theoretical predictions for the form factors, we determine the magnitude of the CKM 
matrix element \Vub\- We find \Vub\ = (3.3 ± 0.4stat. i 0.2syst.lo;4FF) x 10~^, where the last error is 
due to normalization of the form factor. 
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1 INTRODUCTION 



The success of the B Factories has significantly improved our knowledge of the CP violation in 
the quark sector. In particular, the angle /? of the Unitarity Triangle has been measured to a 
5% accuracy from time-dependent CP asymmetries in 6 — ccs decays pQ. On the other hand, 
experimental determination of the other two angles and of the lengths of the two sides (with the 
third side normalized to unit length) have yet to achieve comparable precision. One of the two sides, 
the one opposite to the angle /3, is of particular interest. The uncertainty of this side is dominated 
by the smallest element \Vub\, which is known to about 11% precision. Improved determination of 
\Vub\ therefore translates directly to a more stringent test of the Standard Model. 

Charmless semileptonic decays of the B mesons provide the best probe for \Vub\- Measurements 
can be done either exclusively or inclusively, i.e., with or without specifying the hadronic final state. 
Since both approaches suffer from significant theoretical uncertainties, it is important to pursue 
both types of measurements and test their consistency. 

The exclusive B X^ii^ decay rates are related to \Vub\ through form factors. In the simplest 
case of -B — > irii/, the differential decay rate (assuming massless leptons) is given by 



|2 



a,^ (1) 

where is the invariant-mass squared of the lepton-neutrino system, and is the momentum 
of the pion in the rest frame of the B meson. The form factor (FF) f+{q^) is calculated with a 
variety of theoretical models. In this paper, we consider recent calculations by Ball and Zwicky [2] 
based on light-cone sum rules (LCSR) and by the HPQCD (2j and FNAL |4j Collaborations based 
on unquenched lattice QCD (LQCD). The LCSR and LQCD calculations provide the form factor 
with reliable uncertainties only in limited ranges of . It is therefore necessary to extrapolate the 
calculated form factor using empirical functions, or to measure the partial decay rates Ar(i3 — > it(.v) 
with appropriate cuts on , typically chosen as q^ < 16GeV^ and q"^ > 16GeV^ for use with the 
LCSR and LQCD calculations, respectively. 

Measurements of the partial branching fractions AB{B nil/) have been reported by CLEO [H], 
Belle 0, and by BABAR Y;. The CLEO and BABAR measurements reconstructed B — > Triu events 
by inferring the neutrino momentum from the missing momentum; the Belle measurement used 
B mesons recoiling against another B meson reconstructed in semileptonic decays. BABAR has 
also reported a measurement of the total branching fraction B{B — > nil') using the recoil of fully- 
reconstructed hadronic B decays 0. 

In this paper, we report preliminary results from a study of the B^ — > n~i^u decay, using an 
event sample tagged by B^ — > D^*^^i^l7 decays.^ A similar study of the B^ — > n^i'^i' decay is 
reported in a separate paper [Uj. 



2 THE BABAR DETECTOR AND DATA SAMPLES 

This measurement uses the e^e~ colliding-beam data collected with the BABAR detector JOI at 
the PEP-II storage ring. Charged particles are measured by a combination of five-layer silicon 
microstrip tracker and a 40-layer central drift chamber, both operating in a 1.5 T magnetic field. 
A detector of internally refiected Cherenkov light provides charged kaon identification. A CsI(Tl) 

^Charge-conjugate modes are implied throughout this paper. 
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electromagnetic calorimeter provides photon detection and, combined with the tracking detectors, 
electron identification. The instrumented flux return of the magnet identifies muons by their 
penetration through the iron absorber. 

The data sample analyzed contains 232 million e^e~ — > BB events, where BB stands for 
B^B~ or B^B^. It corresponds to an integrated luminosity of 211 fb~^ on the T(4S') resonance. In 
addition, a smaller sample (22fb~^) of off-resonance data recorded at approximately 40MeV below 
the resonance is used for background subtraction and validation purposes. 

We also use several samples of simulated e"^e^ BB events for evaluating the signal and 
background efficiencies. Charmless semileptonic decays B — > XuH-i' are simulated as a mixture of 
exclusive channels (X„ = vr, r/, r/', />, and uj) based on the ISGW2 model and non-resonant 
B — > Xu^v decays jl2| with hadronic masses above 2771,^. For the signal channels, we give weights 
to the simulated events in order to reproduce the distribution predicted by the recent LCSR 
and unquenched LQCD 1^ calculations as well as the ISGW2 model. 

3 ANALYSIS METHOD 

The analysis method we use for event selection and signal yield extraction has been developed blind, 
i.e., without using the signal sample in the data. The procedure described in the following sections 
has been chosen and optimized using MC simulation to obtain the largest expected statistical 
significance of the partial signal yields in the three bins defined in Section [3.21 

The outline of the analysis is as follows: We look for combinations of a or D*~^ meson and 
a lepton (e~ or that are kinematically consistent with B^ D^*^~^i~17 decays. For each such 
B candidate, we define the recoil side as the tracks and calorimeter clusters that are not associated 
with the candidate. We search in the recoil side for a signature of a B^ it~£~^v decay. We take 
advantage of the simple kinematics of the B^ ir^i^v process to define discriminating variables, 
and extract the signal yield from their distributions in three bins of q"^. Finally we calculate the 
total and the partial branching fractions using the signal efficiencies predicted by a Monte Carlo 
(MC) simulation. We correct for the data-MC efficiency differences using a control sample in which 
both B mesons decay to tagging modes. 

3.1 Event Selection 

We search for candidate BB events in which one B meson decayed as B^ — > D^*')~^i~T', where the 
D*~^ meson is reconstructed in the D*~^ — > D^n^ and D^n^ channels. The D mesons are recon- 
structed in the —>■ K~it~^ , K^tt^tt^tt^ , K^tt^t:^, i^gVr+Tr", and D'^ K^tt^tt^ channels. 
The widths of the signal regions around the nominal D-meson masses are between ±15MeV and 
=b30MeV, which correspond approximately to =b3o" of the mass resolution. We also define sideband 
regions, evenly split below and above the signal regions, which are used to subtract the combina- 
torial background. The sideband regions are chosen to be 1.5 times wider than the signal regions. 
The difference between D*^ and D masses must be within ±3 MeV of the nominal values. 

The candidates are combined with an identified electron or muon to form a V 

candidate. The lepton must have a center-of-mass momentum^ p| > 0.8 GeV. If the D meson was 
reconstructed with a charged kaon, the kaon charge and the lepton charge must have the same sign. 

For each W D^*'>+l'-u candidate, we remove from the event the tracks and the neutral 
clusters that make up the D^*^^£~V tag. We then search for a B'^ — s- Tr^i'^v candidate in the 

^Variables denoted with a star (a;*) are measured in the T{4S) rest frame; others are in the laboratory frame 
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remaining part of the event. We require an identified lepton with p| > 0.8 GeV, accompanied by 
an oppositely charged track that is not identified as either a lepton or a kaon. To allow for the 
B^-B^ mixing, we do not require the two leptons in a candidate event to be oppositely charged. 

The signal events we try to identify contain two neutrinos, one from each B decay. Four- 
momentum conservation and the invariant masses of the B mesons and neutrinos provide just the 
sufficient number of constraints (8) to determine the event kinematics. Referring to the D'^*^~^1~ 
system as the "F" system, we first calculate the cosine of 6*^^, the angle between and py, as 

• 

The energy E*^ and momentum p*^ of the B meson are known from the beam energies and the B^ 
mass. Equation ((JJ assumes that a B^ D^*^^£~T7 decay has been correctly reconstructed, and 
the only undetected particle in the final state is the neutrino. If that is the case, cos^sy should be 
between —1 and +1 within experimental resolution. If the tag has been incorrectly reconstructed, 
Equation (j2j does not give a cosine of a physical angle, and cos 6 by is distributed more broadly. 
Analogously, we can calculate the cosine of 9Bn£, the angle between p^ and p*^, as 

cos Ub-kI - TTT—* ■ 

This variable, again, should be between —1 and +1 for the signal events, and distributed broadly 
for the background. 

The momenta of the two B mesons must be back-to-back in the center-of-mass frame. Given 
cos 9 BY, cos 9 B-n-e, and the directions of py and p*^, we can determine the direction of the B 
momenta up to a two- fold ambiguity. Denoting (pB to he the angle between p^ and the plane 
defined by py and p*^, we find 

2 COs"^ 9 BY + COs'^ 9 BttI + cos 9 BY cos 9 Bni cos J , , 

cos (pB = , (4) 

sm 7 

where 7 is the angle between py and p*^, as shown in Figure ^ It should be noted that cos^^s 
satisfies 

cos'^ (t>B > cos^^sy and cos^(/>b > cos^ 9 b-kI (5) 




Figure 1: Angles (jjB-, 9by, 9btt£, and 7. The B and B momenta (dashed red arrows) are out of the 
plane defined by the Y and ttI momenta (black arrows). 
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by construction. Equation @ assumes that cos^^y and cosOb-kI are correct, i.e., both the tag-side 
and the signal-side of the event have been correctly reconstructed with only one neutrino missing 
on each side. In that case, the variable cos^0b must be between and 1. If any part of the 
reconstruction is incorrect, cos^ 4>b is no longer a cosine squared of any physical angle, and its 
distribution spreads beyond 

We use cos'^ as the principal discriminating variable for this analysis. Since cos^ <j)B uses the 
information of the complete event and is strongly correlated with cos 9 by and cos 6 b-kI, we apply 
only very loose cuts, | cos^sy | < 5 and | cos 6*^77^1 < 5, in the event selection. 

If a signal event has been fully reconstructed, no other particles should be present. In reality, 
such an event often contains extra photons, some of which come from vr'^s and/or photons from 
decays of D* or heavier charmed mesons. Although we use the D*^ — > D^i:^ decay in the tags, the 
efficiency for reconstructing the soft tt" is low. We do not use the D*^ — > D'^^ decay due to a poor 
signal-to-background ratio. We identify the photons that may have come from D* — > Dir^ and 
decays by combining them with the D meson candidate; if the combination satisfies m£)^ — mo < 
150 MeV and cos^^y < 1.1, where Y' stands for the Dj£ system, the photon is considered as a 
part of the D^*^£u system and is removed from the recoil system. Another source of extra photons 
is Bremsstrahlung in the detector material. If either lepton is an electron, we identify and remove 
the Bremsstrahlung photons based on their proximity to the direction of the electron track. 

At this point, we require that the event contains no charged tracks besides the D^*^^£~V 
and tt~£'^i' candidates. We further require that there be no residual photons with laboratory- 

frame energy E^^ > 100 MeV. 

A few additional cuts are applied to improve the signal-to-background ratio by rejecting specific 
types of known background events. In order to suppress non-BB background, we require that the 
ratio i?2 of the second and zeroth Fox- Wolfram moments computed using all charged tracks 
and neutral clusters in the event, be smaller than 0.5. The invariant mass of the two leptons, if they 
are oppositely charged, must not satisfy 2.95 < r^g+g- < 3.15 GeV or 3.05 < m^+^- < 3.15 GeV 
in order to suppress the background due to J/^p £~^i~ decays. We also calculate the invariant 
mass of the 'ir~^i~ system assuming that the pion was a lepton of the same species as the identified 
lepton, and require this mass to be outside 3.06-3.12 GeV. 

After all the selection cuts, a small fraction of events contain more than one candidate; simulated 
signal events that pass the selection contain 1.14 candidates on average. When there are multiple 
candidates in an event, we select the candidate with the smallest value of | cos 6by\- The candidates 
in the D-mass sidebands are included in the selection procedure. 

3.2 Signal Yields 

We find in the on-resonance data 966 it 31 and 725 it 34 candidate events in the D-mass signal 
region before and after the sideband subtraction, respectively. Table ^ summarizes the numbers of 
candidate events for each decay channel of D^*^^ . 

We extract the signal from the cos^ cf)B distribution in three (f' bins, namely, (f' < 8GeV^, 
8 < < 16GeV^, and > 16GeV^, as shown in Figure |2 For each event passing the selection, 
we calculate q'^ by 

= {Pl+Vuf ~ ipB -PTrf, (6) 

where p_B is the approximate B four- momentum defined, in the center-of-mass frame, as 

p*^ = {E*B,p*B)^(^^^^,0,0,Oy (7) 
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Table 1: Numbers of candidate events that enter the signal-yield fit from the on-resonance data. 
The D-mass sidebands have been subtracted. 



decay channel 


Events 


U yj\ 7T 7T j 


OOU.O It ZO.Z 




89.0 ± 9.9 




102.3 ±11.7 




172.7 ±14.8 




17.7 ± 5.2 




7.3 ± 3.4 


Total 


725.3 ±33.6 



In other words, the center-of-mass motion of the B meson is ignored. Since the B momentum in 
the center-of-mass frame is small, the impact of this approximation is small, as will be shown in 
Section f3. 31 It is worth noting that the experimental input to Equation ® is the pion momentum 
and not the lepton momentum. As a result, electron energy loss due to unrecovered Bremsstrahlung 
has no impact on the resolution. 

The backgrounds in this measurement are handled in three groups. First, the combinatoric 
background for the D mesons is subtracted using the D-mass sidebands. The remaining back- 
grounds are mostly BB events, and are separated from the signal using the cos'^ (ps distribution. 
The possible contribution from the non-BB background is estimated, based on the cos^ (pB distri- 
bution of the off-resonance data events that pass the event selection, to be smaller than 1.0 event 
of the total yield, and is included in the systematic error. 

The raw signal yield is extracted in each q'^ bin by a simple binned fit of the cos^ (pB distribu- 
tion of the on-resonance data to the weighted sum of the signal and BB background distributions 
from the MC simulation. The sources of the BB background can be -B — > Xulu decays, B — > Xcii' 
decays, and other (hadronic) B decays. The cos^ (pB distribution for the signal events peaks between 
and 1, while that of the background is broad with a gradual fall off toward large values of cos^ (pB- 
The fall-off is faster for smaller g^, as it can be seen in Figure EJ Since the cos^ cpB distributions 
for the BB background from various sources are quite similar, we fix their relative abundances in 
the fit, and later vary them within their systematic uncertainties. The fit therefore has two free 
parameters: the normalization of the signal and the normalization of the BB background. 

To maximize the statistical sensitivity, the first bin of the cos^ (pB histogram should be at least 
as narrow as the signal peak. At the same time, each bin should contain sufficient entries to prevent 
the fit from becoming biased. We use variable bin sizes to satisfy the two requirements. The 
bin boundaries are chosen to be cos^ </>b = 0, 1, 2, 4, 7, and 12. 

TablelSlsummarizes the signal yield obtained from the fit of the on-resonance data, the values 
(for 3 degrees of freedom) and the corresponding probabilities. The combined probability of all 
three bins is 0.033. The poor x^ value of the last q'^ bin is caused by the large number of events 
in 1 < cos^ (pB < 2. Moving the first bin boundary from 1.0 to 1.5 improves the x^ and increases 
the yield significantly. We can find no reason for this behavior other than statistical fluctuation. 
We therefore retain the result obtained with the original binning, which was chosen before the data 
were unblinded. 
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Table 2: Partial signal yields, values and probabilities in bins of . The errors are statistical. 



bin 


Yield (events) 




Prob. 


g2 < 8 GeV^ 


26.3 ±8.7 


4.6 


0.20 


8 < g2 < 16 GeV^ 


21.2 ±9.2 


1.9 


0.59 


> 16GeV2 


14.2 ±8.8 


11.7 


0.008 
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Figure 2: The cos2,/.b fit of the data for q"^ < 8GeY^ (left), 8 < < leGeV^ (middle), and 
q^ > 16GeV^ (right). The rriD sidebands have been subtracted. The points with error bars are the 
on-resonance data. The histograms are signal (yellow) and BB background (cyan). 



3.3 Signal Efficiencies 

In order to derive the partial branching fractions from the signal yields, we need the signal efficiency 
in each q'^ bin. Since the measured value of from Equation © in principle differs from the true 
value, we define the efficiency Eij as the probability, averaged over electrons and muons, of a 
B^ Tr~i^iy event whose true value belongs to the j-th bin to be found in the i-th measured-g^ 
bin. With this definition, the signal yield obtained in the i-th reconstructed-g^ bin is expressed as 

N, = 2j2eijABjNB, (8) 

j 

where the factor of two comes from using both electrons and muons, ABj is the partial branching 
fraction in the j-th q^ bin, and Nb is the number of B^ mesons in the data sample. Using the 
r(45) B^B^ branching fraction /oo = 0.488 ±0.013 the number Nb equals 2fooN^^, where 
is the number of BB events in the data sample, and the factor of 2 comes from having two 
B mesons in each event. 

We use the Monte Carlo simulation to estimate Eij and correct for the known data-MC differ- 
ences as 

pdata 

N^ = '^-^T.^^''^^^Nb. (9) 

j 

Here we assumed that a single data-MC efficiency ratio gdatay^MC ^^^^ applied to all q"^ bins. 
This is reasonable so long as the ratio is close to unity, which is in fact the case as it will be shown 
below. 
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The efficiency determined from the MC simulation is 



MC 
ij 



^ 1.142 ±0.065 0.050 ±0.011 0.004 ± 0.004 \ 

0.074 ±0.015 1.232 ±0.071 0.035 ± 0.015 
^ 0.007 ±0.008 0.063 ±0.016 1.350 ± 0.097 y 



10"^ (10) 



The errors are due to Monte Carlo statistics. The selection efficiency averaged over the three q'^ 
bins is 1.32 x 10^'^. In addition to the experimental resolution, final-state radiation (FSR) changes 
the q'^ distribution. All MC samples used in this analysis are generated with PHOTOS jl5j to 
simulate FSR. The bin-to-bin migration due to FSR is predicted by the MC simulation to be less 
than 1.2%, and the full size of this effect is included in the systematic error. 

We evaluate the data-MC difference of the B'^ D^*^^i~V selection efficiencies using the 
double-tag events, in which both B mesons decay to D^*^^£i'. Properties of the D^*^£u tags such 
as the composition of the D^*^ decay channels are similar for the tagged-signal and double-tag 
events. The number of double-tag events is proportional to the square of the tagging efficiency 
after subtracting the small contribution from background. 

The selection criteria for the double-tag events follow the main analysis as closely as possible. 
In each event, we look for two D^*^lv tags that do not share any particles. We remove all particles 
that are used in the two tags and require that there be no charged tracks and no neutral clusters 
remaining in the event. 

After subtracting the Z?-mass sidebands, we find 1073.4 double-tag events in the on-resonance 
data. Figure 131 shows the cos^ 4)b distribution of the selected events. The 'signal' in this case consists 
of the B^B^ events in which the two B mesons decay into D^*^iv and are correctly reconstructed 
as two tags. A small fraction of B^B^ events with two D^*^£v decays are incorrectly tagged, i.e., 
wrong combinations of particles are selected as the two tags. Other sources of background events 
include B ^ t ^ i cascade decays and lepton misidentification. 



Preliminary 




Figure 3: Distribution of cos^ (pB for the double-tag events in the data (error bars) and in the MC 
simulation (histograms). The niD sidebands have been subtracted. Colors indicate events with two 
correct tags (yellow), BB events with two B D^*^lv decays that were not correctly reconstructed 
(red), and other background (cyan). 
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We extract from the study of the double-tag events the efficiency correction factor 

g-data 

= 1.000 ± 0.047. 

The stated error includes both the statistical and systematic uncertainties. For the latter, we 
considered a) the difference in the results when the selection criteria are relaxed to allow presence 
of neutral clusters, b) the residual background after the D-mass sideband subtraction, due to 
possible non-linearities in the backgrounds vs. D mass, and c) the uncertainties in the exclusive 
B — > Xcii' branching fractions. The relaxed criteria used in a) increase the double-tag yield and 
the background by approximately 50% and 100%, respectively. 

4 SYSTEMATIC UNCERTAINTIES 

The significant sources of systematic uncertainties and their impact on the measured total and 
partial branching fractions are summarized in Table |21 

For the B — > nil' form factor, we use the Ball-Zwicky calculation for our central values, and con- 
sider the differences between that and the two LQCD calculations as the systematic uncertainties. 
The form factor affects the branching fractions through the dependence of the signal efficiency. 



Table 3: Fractional systematic errors on the total and partial branching fractions. 







0"AB 


/AB (%) in g2 


bins 


Systematics 


ats/B (%) 


< 8 GeV^ 


8-16 GeV^ 


> 16 GeV^ 


B -Kiv FF 


-h0.2/-0.0 


±0.2/-0.7 


±0.7/-0.3 


±1.6/-0.3 


B plv FF 


-F0.5/-0.0 


±0.0/-0.2 


±0.5/-0.0 


±2.6/-0.0 


B{B plv) 


-F2.9/-3.4 


±0.1 


±3.7 


±7.8/-10.2 


B{B XJu) 


±2.8 


±4.8 


±9.9 


±13.0 


B{B XJu) 


±1.7 


±0.1 


±1.9 


±5.6 


Shape Function 


±0.8/-l.l 


±0.1 


±0.6 


±2.8/-4.1 


D(*ku Tagging 


±4.7 


±4.7 


±4.7 


±4.7 


Tracking 


±1.6 


±1.6 


±1.6 


±1.6 


Electron ID 


±1.1 


±1.2 


±1.1 


±0.9 


Muon ID 


±1.4 


±1.2 


±1.4 


±1.6 


resolution 


NA 


±2.0 


+0.6/-0.8 


±2.9 


Final-state radiation 


NA 


±1.2 


±1.2 


±1.2 


cos^ (j)B for backgrounds 


±8.8/-9.1 


±5.6/-1.6 


+10.1/-14.2 


+14.0/-17.9 


l^on-BB background 


±1.6 


±1.6 


±1.6 


±1.6 


Combinatoric D 


±0.8 


±0.8 


±0.8 


±0.8 




±1.1 


±1.1 


±1.1 


±1.1 


/oo 


±2.7 


±2.7 


±2.7 


±2.7 


MC statistics 


±3.5 


±5.8 


±6.1 


±7.6 


Total 


±12.2/-12.5 


±11.5/-10.2 


+17.2/-19.9 


+24.1/-27.4 
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As a result, only the shape and not the normalization of the form factor is relevant at this stage, 
while the normalization becomes important in the determination of \Vub\ as discussed in Section 1^1 
The B p£u decays are significant sources of background at large q'^. We vary the branching 
fractions as ^ p-l+u) = (2.69l[J:^^) x 10"^ and B{B+ p^l+u) = (1.45l[J:^?) x 10-^ based 

on the measurements |16| and isospin symmetry. We also compare results using the B — > piv form 
factors calculated by Ball and Zwicky _2^, Melikhov and Stech ^j, and UKQCD |18j . 

The branching fractions for the B — > X^H-y and B Xulu decays also significantly affect the 
backgrounds. We use the latest measurements for the branching fractions. Where appropriate, 
we combine the B^ and B^ branching fractions assuming isospin symmetry. The shape function 
parameters used for the simulation of the non-resonant B — > Xuif decay are varied according to 
Ref. Uni. 

The efficiency for the D^*^ii' tagging has been discussed in Section 13.31 The systematic un- 
certainties in the track reconstruction and lepton identification efficiencies have been derived from 
studies of independent control samples. We vary the amount of migration between the bins 
due to resolution, given by the off-diagonal components in Equation (fTn|) . by ±50%. We assign an 
additional ±1.2% error on the partial branching fractions for the g^-bin migration due to final-state 
radiation, which is simulated using PHOTOS |l2j • 

The largest source of systematic error is the shape of the cos^</>b distribution for the BB 
background events. We studied it using several control samples that are depleted of the irii' 
signal, obtained, for example, by requiring one extra charged track or neutral cluster remaining in 
the event. The cos^ (f)B distributions of the control samples agree between the data and the MC 
simulation within the available statistics. We assign the systematic errors based on the statistical 
uncertainties of these tests. 

We used the off-resonance data sample to set an upper limit on the residual non-BB background. 
MC simulation was used to determine the size of the combinatoric background that remains after 
the D-mass sideband subtraction, due to the non-linearities in the backgrounds vs. reconstructed 
D mass. 

The number of BB events in the on-resonance data sample is known to ±1.1%. We use /oo = 
0.488 ± 0.013 ,14^ as the T{AS) B^B° branching fraction. Limited statistics of the MC samples 
affects the measurement primarily through the estimation of the signal efficiency e}j^ given in 
Equation ((TU)) . 

In addition to the studies discussed above, we perform a large number of crosschecks and tests 
of cut-value dependences. We investigate all cases in which the variations exceed the expected 
statistical fluctuations, and find no indications of systematic problems. 

5 RESULTS 

From the signal yields and the efficiencies evaluated in Section IHl we extract the following prelimi- 
nary result for the total branching fraction: 

B{B^ ^ TT-i+u) = (1.03 ± 0.25stat. ± 0.13syst.) x 10"^ 

The partial branching fractions and their statistical errors are given in Table [l] Note that the 
errors in the partial branching fractions are negatively correlated because of the small migration 
across bins, so adding them in quadrature does not give the error in the total branching fraction. 
We extracted the branching fractions for each of four signal form-factor calculations (Ball-Zwicky, 
HPQCD, FNAL, ISGW2); the differences are smah. 
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Table 4: Preliminary results for the partial and total tt i^u branching fractions. The errors 

are statistical. 







10^ 




Signal FF 


< SGeV^ 


8-16 GeV^ 


> 16GeV2 


B X 10^ 


Ball-Zwicky [2] 


0.483 ±0.166 


0.337 ±0.162 


0.209 ±0.140 


1.029 ±0.253 


HPQCD 


0.480 ±0.165 


0.339 ±0.163 


0.211 ±0.142 


1.031 ±0.254 


FNAL m 


0.483 ±0.166 


0.337 ±0.162 


0.209 ±0.141 


1.029 ±0.253 


ISGW2 11 


0.477 ±0.165 


0.339 ±0.163 


0.215 ±0.147 


1.032 ±0.256 



xlO' 



V 0.6F 
t 0.5F 



<i 



0.4 1 
0.3 j 
0.2; 
0.1; 

0- 



Ball-Zwicky 

— HPQCD '04 
— ■ FNAL '04 

— ISGW2 



BaBar 

Preliminary 



q2<8 



8 < q < 16 q2 > 16 (GeV ) 



Figure 4: The partial branching fractions AB{B^ — > vr i^u). The points with error bars are the 
measurements; the thick parts of the error bars indicate the statistical errors. The histograms are 
predictions by Ball-Zwicky ^ (red), HPQCD (green), FNAL m (dashed blue), and ISGW2 [III 
(black) calculations, each scaled to the measured total branching fraction. Thus only the de- 
pendences are relevant in this comparison. The Ball-Zwicky (red) and FNAL (dashed blue) lines 
overlap with each other. 



Figure 13 compares the measured partial branching fractions with the dependence predicted 
by Ball-Zwicky, HPQCD, FNAL, and ISGW2 calculations. The calculations are normalized to 
the measured total branching fraction. Although the measured AB values depend slightly on the 
form-factor calculation, the differences are too small (< 2.5%) to be noticeable on this plot. 

Table El summarizes the measurements of B{B — > Triv) by the BABAR collaboration. Assuming 
isospin symmetry and the ratio of B lifetimes tq+/tqo = 1.081 ± 0.015 J^, the measurements 
agree with each other with = 10.3 for 5 degrees of freedom, which corresponds to a one-sided 
probability of 7%. 

Prom the measurement of the partial branching fractions AB described in this paper, we extract 
\Vub\ using 
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Table 5: BABAR measurements of B{B — > Triu). All results are preliminary. The last row shows 
this measurement and the B~^ — > ir^l^u measurement in Ref. fS]. 



Technique 


B{B^ - 




TT-i+l^) X 10^ 


B{B+ 




7r0£^ 


'I' 


) X 10^ 


Neutrino reco. 


1.41 


± 


0.17 ±0.20 


0.70 


± 


0.10 


± 


0.10 


Hadronic tag |S1 


0.89 




0.34 ±0.12 


0.91 


± 


0.28 


± 


0.14 


Semileptonic tag 


1.03 




0.25 ±0.13 


1.80 


± 


0.37 


± 


0.23 



Table 6: Preliminary results of \Vub\ extracted from the measured partial (first three rows) and 
total (last three rows) branching fractions and form- factor calculations. 



FF calculation 


Q 


^ range 


AC(ps-i) 




\Vub\ 


(10-3) 




Ball-Zwicky |21 


< 


IGGeV^ 


5.44 


± 


1.43 


3.1 


± 0.4stat. 


± 


0-2syst. 


-fO.5 

-0.3FF 


HPQCD 


> 


leCeV^ 


1.29 


± 


0.32 


3.3 


± 1-lstat. 


± 


0.5syst. 


+0.5 

-0.3FF 


FNAL U 


> 


IGGeV^ 


1.83 


± 


0.50 


2.7 


± 0.9stat. 


± 


0.4syst. 


+0.5 
-O.3FF 


Ball-Zwicky (2j 




full 


7.74 


± 


2.32 


2.9 


± 0.4stat. 


± 


0-2syst. 


+0.6 
-O.4FF 


HPQCD 




full 


5.70 


± 


1.71 


3.4 


± 0.4stat. 


± 


0-2syst. 


+0.7 
-O.4FF 


FNAL 




full 


6.24 


± 


2.12 


3.3 


± 0.4stat. 


± 


0.2syst. 


+0.8 
-O.4FF 



where r^o = 1.536 ± 0.014 ps ^19^ is the B lifetime, and A(" is defined as 

^C = ^/, l/+(.^)lW- (12) 

^min 

To minimize the theoretical error on \Vub\, the range of q'^ should correspond to the region in which 
the form-factor calculation is most reliable: < 16 GeV^ for LCSR and q"^ > 16 GeV^ for LQCD. 
In order to extract \Vub\ from the total branching fraction B, the form factor must be extrapolated 
to the full range of q^. This is done in Kefs. IS il using empirical functions, and additional 
uncertainties are assigned for the extrapolation. 

Table ini summarizes the values of \Vub\ extracted from the measured partial and total branching 
fractions. The last errors on \Vub\ come from the uncertainties in A^, which in turn come from 
the uncertainties of the normalization of the form-factor calculations. The precision of the results 
obtained using the LQCD calculations |31I5 9^ > 16 GeV^ are limited by the large statistical error 
of the measured partial branching fraction. Using the total branching fraction reduces experimental 
errors on \Vub\ at the cost of increased theoretical uncertainties due to the extrapolation of the form 
factor. 

Instead of averaging results based on different theoretical calculations, we report the value of 
\Vub\ obtained from the total branching fraction based on one of the LQCD calculations jl], 

\Vub\ = (3.3 ± 0.4stat. ± 0.2syst.lo:4FF) X lO^^. 

as a representative result. This result lies between the results based on the other two calculations, 
and carries the most conservative theoretical uncertainty. 



18 



6 SUMMARY 



Using event samples tagged by D^*^^£ u decays, we obtain the exclusive branching fraction 

B{B^ 7r~£~^i'). The preliminary result for the total branching fraction is 

13{B° ^ 7r~e+u) = (1.03 ± 0.25stat. ± O.lSsyst.) x 10~^ 

We also obtain the partial branching fractions in three bins of q'^ of the lepton-neutrino pair. The 
preliminary results are 

[ (0.48 ± 0.17stat.lo:o5syst.) X lO'^ < 8GeV2, 
Ae(S° ^ TT-i+i^) = I (0.34 ± 0.16stat.l[l:[!?syst.) X 10"^ 8 < q^ < 16GeV2, 
[ (0.21 iO.Mstat.laSesyst.) x 10"^ q^ > IGGeV^. 

Using the measured total branching fraction and a form- factor calculation based on lattice QCD [1], 
we extract 

\Vub\ = (3.3 ± 0.4tat. ± 0.2syst.loiFF) X lO^^. 

where the last error is due to normalization of the form-factor. We also use other recent calculations 
of the form factor [2113] and find values of \Vub\, given in Table El that are consistent within the 
experimental and theoretical uncertainties. 
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